The rigidity transition in polymer melts with van der Waals interactions 
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We study the onset of rigidity near the glass transition (GT) in a short-chain polymer melt mod- 
elled by a bead-spring model, where all beads interact with Lennard- Jones potentials. The properties 
of the system are examined above and below the GT. In order to minimize high cooling-rate effects 
and computational times, equilibrium configurations are reached via isothermal compression. We 
monitor quantities such as the heat capacity Cp, the short-time diffusion constants T>, the viscos- 
ity rj, and the shear modulus; the time-dependent shear modulus G{t) is compared with the shear 
modulus jj, obtained from an externally applied instantaneous shear. We give a detailed analysis of 
the effects of such shearing on the system, both locally and globally. It is found that the polymeric 
glass only displays long-time rigid behavior below a temperature Ti, where Ti < To- Furthermore, 
the linear and non-linear relaxation regimes under applied shear are discussed. 



PACS numbers: 61.25.Hq, 66.20.-|-d, 64.70.Pf, 65.60.4-a 



I. INTRODUCTION 
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The onset of rigidity in disordered systems has been the 
focus of a number of studies over the last decades 0, |3| ■ 
At the rigidity transition, a resistance to shear emerges, 
characterized by a non-zero shear modulus (above the 
transition), and a diverging viscosity (below the tran- 
sition). In an important class of systems the onset 
of rigidity can be understood in terms of percolation 
[IIIHIEIEIIHH' Connectivity percolation coin- 
cides with the onset of an entropic component to rigidity 
0, IMBji of the type traditionally associated with rub- 
ber [13. Mechanical rigidity requires a rigid percolat- 
ing backbone. If the interaction forces are central this 
means a multiply connected backbone, and consequently 
rigidity would appear at a higher concentration of bonds 
than connectivity percolation j^J. This physical picture 
provides a well-defined approach to understanding the 
behavior of the viscosity and the shear modulus in the 
neighborhood of the onset of entropic rigidity, and has 
been the focus of several recent studies 0, H, IE S IS 13 ■ 
The onset of the mechanical rigidity arising from covalent 
bonds has been extensively studied in network glasses at 
temperatures well below the glass transition (GT). In a 
mean-field argument, this component of rigidity sets in 
when the density of floppy modes approaches zero, or the 
number of constraints exceeds the number of degrees of 
freedom. For bond-bending networks, this means a mean 
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In this paper we consider the onset of rigidity in a poly- 
mer melt which occurs as the temperature is lowered be- 
low the GT. There is no obvious length scale emerging 
a priori as the chains bind together under the attractive 
van der Waals interactions between monomers forming 
disordered structures. As we shall see, the region around 
the glass transition, where the viscosity is expected to 
diverge and a finite shear modulus to emerge, becomes 
more complex than in percolation driven systems. Where 
precisely the onset of rigidity occurs is an open question 
[l5| . There is a body of literature on the behavior of 
the viscosity near the glass transition which will be used 
to discuss our results. There is however much less infor- 
mation about the onset of rigidity as measured by the 
emergence of a shear modulus. 



The article proceeds as follows. Section II provides 
some general, conceptual background information re- 
garding the GT, and in particular, the behavior of the 
viscosity. Section III discusses our model and the compu- 
tational tools that we have used. Section IV presents the 
main results obtained from these simulations and some 
brief explanations. In addition to the viscosity and the 
shear modulus we have calculated the diffusion constant 
and the heat capacity. Section V, before the Conclu- 
sion, discusses the results within the current framework 
of polymeric glasses and provides some insight into the 
theoretical concepts discussed. 
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II. RIGIDITY AND THE GLASS TRANSITION 

For a polymer melt, in contrast to the rigidity that 
sets in with increasing crosslink density, lowering the 
temperature drives the system reversibly towards a ho- 
mogeneous glass phase with no obvious diverging length 
scale. Attempts have been made to invoke the percola- 
tion of the slow regions near the glass transition [l7l Il8l | . 
And also in a particular model where rigid bonds can be 
formed upon cooling, a link is made to a mechanical rigid- 
ity transition based on constraint counting j3| • We 
will focus on the behavior of the viscosity and the shear 
modulus and analyze their critical behavior in a simple 
numerical model of a polymer melt made of short freely 
jointed chains. Discussions will be made using existing 
models of the physical properties of the system near the 
glass transition. 

Historically viscosity has played an important role in 
analyzing glasses. The glass transition was associated 
for practical purposes with a specific number for the vis- 
cosity: 77 = 10"'^^ poise or 10"'^^ Pa-s, where the material 
could be considered not to flow during experimental time 
scales. 

The glass transition is known as a "pseudo second- 
order transition" , since the discontinuities are not sharp 
and occur over a range of temperatures As well, a 
glassy system is not actually in equilibrium, because ag- 
ing phenomena and flow have a non-negligible effect at 
long times. At best, we refer to a metastable equilibrium 
which, on our timescales (as in the case of short experi- 
mental times), can be treated as an equilibrium system. 

A commonly-used approach to the glass transition has 
been mode-coupling theory (MCT), derived from micro- 
scopic density fluctuations within a liquid A 
dynamic phase transition is predicted at approx. 1.1-1.3 
times Tg, involving a divergence in the zero shear-rate 
viscosity rj, 



r]={T - TmcY 



(1) 



which could be masked by thermal activation effects [2J| . 
More importantly, idealized MCT (which does not take 
into account "hopping" processes, for instance) predicts 
an ergodicity-breaking at T = Tmc [13 ■ It has also been 
proposed that there is, in fact, some type of "crossover 
point" to the Vogel-Fulcher-Tamman (VFT) regime (see 
[2z| and Refs. therein). The latter is a more success- 
ful framework for analyzing glass-forming liquids and is 
defined by the VFT law, an empirical equation of the 
form 



77 = 77oo exp 



T-To 



(2) 



The VFT temperature Tg is expected and has been ob- 
served to be very close to the Kauzmann temperature 
discussed in Refs. [l^H^. Tk is the temperature which 



defines a hypothetical state with zero configurational en- 
tropy or the formation of an "ideal" glass. Such a state 
is never reached because of the large slowing-down of the 
system at Tg- Impressive fits to the VFT equation have 
been achieved ^26^ 2Tj above the glass transition. These 
fits also provide a means to characterize the fragility of 
glasses: the closer Tg is to Tq the more fragile the glass 
j28, 29]. Nevertheless, the range of applicability of the 
VFT equation is under contention; some research indi- 
cates that it cannot go as low as Tg, while other studies 
suggest that it is best applied close to the GT (see Ref. 
|25| for a review). 

A dynamical scaling approach was recently proposed 
to explain the divergence in relaxation time below the GT 
[2^. The GT is characterized by dynamical heterogene- 
ity, in that only some of the material is able to move due 
to a lack of free volume. At a critical temperature Tg, 
there is not enough free volume for any particles to move 
and local motion is effectively prohibited. This approach 
assumes percolation-like universal behavior among glass- 
forming liquids with an additional activation energy term 
of the following form 



V 



T-Tg 
Tg 
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(3) 



Unlike the MCT approach, the above equation yields a 
critical temperature for viscosity divergence Tg below 
Tg- While experimental data from polymeric liquids 
is well-described without the activation energy contribu- 
tion, the extra term ensures the same universality class 
for all glass formers [28l |. 

One can also consider the GT in terms of the "en- 
ergy landscape" of the system, a hypersurface punctu- 
ated by local minima of varying depth, where the shape 
depends on the volume and the sampling on the temper- 
ature [2^ [33, 1^ ■ In terms of energy, Tg can be seen 
as a temperature at which the system is "trapped" in a 
relatively deep local minimum and the free energy barri- 
ers effectively prevent further exploration of phase space. 
Polymeric glasses do not behave like other, more sim- 
ple glasses, such as Si02. Simple glasses have a more 
precisely defined energy landscape, whereas polymeric 
glasses have more configurational possibilities, yielding 
more local minima which can be sampled via local rear- 
rangements. As such, polymers tend to create relatively 
"fragile" glasses, characterized by a non-Arrhenius vari- 
ation of the viscosity with temperature. In terms of the 
glass transition temperature, it has been observed that 
Tg rises with chain length or molecular weight until a 
plateau around M ^ 100 p^. 

In general, glasses display a more complex response to 
low-frequency shearing than purely elastic systems. Non- 
linear mechanical responses are common in glass systems 
after ergodicity-breaking '25'| , and the relaxation of a sim- 
ple glass will depend on local rearrangement. In polymer 
systems there is an additional component arising from 
the extension of the chains. This, as well as energy land- 
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scapes, will be useful tools in interpreting the mechanical 
response within our glassy polymer melt. Some aspects 
of relaxation below the GT are well understood [23, El| 
in terms of aging subsequent to a temperature quench 
or in response to a mechanical deformation, after which 
glasses can relax fairly quickly. However this has not, 
to our knowledge, been examined using an instantaneous 
simple shear experiment. 



III. MODEL 



A. The bead-spring model 

We have used the well-known bead-spring model based 
on work by Kremer and Grest [3^ and more recently ap- 
plied with success by the research group of Kurt Binder 
et al. miiUll iHl lli. Their work has focused on 
issues such as a and /3-relaxation, cooling-rate depen- 
dencies and cage effects on approaching Tg from above. 
In this model, the neighboring beads along each chain 
(there are 10 monomers per chain and a total of = 1050 
monomers) interact through the FENE potential 
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(4) 



where Ro — 1.5, k = 30. In addition, all particles interact 
though the truncated Lennard- Jones potential. 
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Tij < 2.5(7 



(5) 

where C is chosen such that the potential is zero at the 
cutoff radius. 

The length is rescaled using a — 1 and the temper- 
ature is expressed in units of eLj/ks- Combining the 
two potentials yields an optimum bond length of 0.96(7, 
which inhibits crystallization by introducing a competing 
length scale. 



B. Approaching the glass phase along an isotherm 

Usually the GT is approached at constant P by low- 
ering the temperature |2^ HE HE ■ In order to 
avoid problems associated with high cooling rates 
and to reduce computational time, we propose a method 
whereby the samples can be "frozen" or not, as the case 
may be, by compression at a desired T. Although we are 
still measuring changes in a given quantity with respect 
to temperature with a constant pressure, we are effec- 
tively using a different thermodynamic path to define the 
state of a sample. There is a v ery complete discussion of 
thermodynamic paths in Ref. [37|. However, it does not 



discuss defining a thermodynamic state using an isother- 
mal path. Our method seems to allow a better explo- 
ration of phase space and consequently, a more physically 
meaningful final state for low temperatures. In a pure 
Lennard-Joncs simulation, for instance, we found that 
crystallization requires less CPU time when compressing 
the system than when simply cooling it j38j |. The tem- 
peratures studied range from 1.2 to 0.2 CLj/ks, allowing 
us to explore the liquid, supercooled and glassy regimes. 
We entered the glass phase both along an isochore and 
an isobar. However, we retained only the simulation re- 
sults along the isobar, as pressure variations along the 
isochore proved somewhat unphysical. The main prob- 
lem was a progressive increase of pressure with lowering 
temperature below Tq (similar results were found in Ref. 
[35|i. This discrepancy effectively amounts to 
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which is clearly forbidden. Along the constant volume 
path, the system goes through what would have been a 
two-phase region in the thermodynamic limit. 

For a given sample, we create an initial "gas" of poly- 
mers at the desired temperature, having ensured their 
separation and random orientations. The chains reach 
their individual equilibrium configurations in an NVT 
ensemble using an algorithm based on the Langevin equa- 
tion, called Brownian Dynamics (BD) in Ref. |39j: 



d Xi 



dxi 



mT- 



,dxi 
~dt 



(7) 



This algorithm allows the system to reach equilibrium in 
our "expanded" volume quickly and realistically, by sim- 
ulating a heat bath in the form of a friction coefficient 
r and by introducing a random force W in the form of 
Gaussian noise. Subsequently, the system is compressed 
to the pressure P — l.QeLj/cr^ to achieve a similar system 
to that of K. Binder's research group pllsl Isl Im ISTj . 
The compression is done in the BD NVT ensemble by 
reducing affinely the system size at a specified rate. The 
largest compression rate of a side of the box, in Lennard- 
Jones units, was around 0.015y^s Lj/m. Using a smaller 
compression rate was shown to cause negligible change 
in quantities such as viscosity and energy. Naturally, 
the rate of compression will inevitably have some effect 
on determining at what temperature Tq the glass be- 
comes "stuck." However, we believe that, by continu- 
ously changing the volume (and thus the precise shape 
of the energy landscape) , the system can eventually find 
a lower energy minimum. The final volume is estab- 
lished at a given temperature for P = 1 using a con- 
stant pressure damped- force algorithm [s^- The system 
is then allowed to evolve in a micro-canonical ensemble 
by MD. We need to run in that ensemble to collect infor- 
mation on the correlators required to calculate viscosity 
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and shear modulus related quantities. The equations of 
motion were integrated with the standard velocity Verlet 

algorithm [s^ and a time step of dt — 0.005 ^ ma'^ /elj- 
During this portion of the simulation, both the pressure 
and the temperature were seen to fluctuate around their 
expected values. 



C. Viscosity and shear modulus from the stress 
autocorrelation function 

The time-dependent shear modulus G{t) measures the 
response of the system to a shear strain e^^ applied at 
t — 0, and is defined as ^3 



Git) 



(8) 



aaf3 is an off-diagonal macroscopic stress tensor element 
and implicitly includes inter and intra-chain interactions. 
G{t) can be calculated from the stress fluctuations in the 
quiescent melt using the fluctuation dissipation theorem 
which states that |l(ll |: 



G{t) 



V_ 

kf 



< crapito)(Tapito + t) P (9) 



In the past, viscosity measurements have been made 
using non-equilibrium molecular dynamics (NEMD) to 
avoid the problem of long time tails of correlation func- 
tions at temperatures approaching Tq '27]. This issue 
was resolved by fitting G{t) to the KWW, or stretched- 
exponential function and integrating to infinity, via Efli | 



G(<) = G(0) cxp 



(13) 



The quality of the fits as well as previous research [23 
suggests that this approach reflects well the behavior of 
the system over long time scales. Although there is a 
certain additional error associated with the fit, it allows 
us to keep a system in its "true", undriven state. The 
errors were large for T very close to Tq {T Tg + 0.02), 
due to the increasing length of the tails. 

The nature of the tail changes significantly below the 
glass transition, and the KWW fit of Eq. is no longer 
appropriate. The slow relaxation of the glass through 
thermally activated pro cesses is better represented by the 
power law form pol |4]| : 



Git) = Ge, + At 



(14) 



where 



In physical terms the power law fit uses the distribution 
of energy barriers explored by the system during the sim- 
ulation time to predict the long time behavior. 
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The sums are over all the particles in the system indexed 
from 1 to TV. G{t) is a very powerful tool: it allows us to 
obtain a zero-shear limit for the shear modulus without 
external deformation, which can have a non-negligible ef- 
fect on the steady state. Non-ergodic glassy materials are 
especially sensitive to configurational changes. Instead, 
G{t) simply uses the microscopic density fluctuations and 
the subsequent response to give us a clear picture of how 
the stress evolves in a system. At < ^ 00, G{t) becomes 
the "equilibrium modulus," Geq and, at t — > 0, Gjt ) be- 
comes Goo , the infinite- frequency modulus jifl, |41| . To 
our knowledge, Geq has never been compared to the usual 
shear modulus, which is measured by a system's response 
to a shear deformation (see Section IIII El below) . The 
zero shear-rate viscosity is obtained by the appropriate 
Green-Kubo formula: 



r^ = j G{t)dt, 


from which one can obtain the relaxation time as 



Goc 



(11) 



D. Short-time diffusion coefficients 

To get an idea of the mobility of the chains and the 
individual particles, we define short time diffusion coef- 
ficients using the Einstein relations, Vm, for individual 
particles i, and T>c'm, for the centers of mass of individual 
chains j: 



'DcM — 



[r. (t) - r, (0)]^ 
6t 



[r, {t) - r, (0)]' 
6t 



(15) 



(16) 



(12) 



Ignoring initial displacements, we simply retain the slope 
of the first approximately linear part of the mean-square 
displacement vs. time graph, calculated around iOOOOdt. 
In the asymptotic limit if the particles are delocalized, 
there would be on a longer timescale a linear regime with 
a smaller slope. If they are localized, (< Ar >)^ would 
be bounded. In the liquid phase, the first linear regime 
essentially reflects the local diffusion of a particle, which 
is limited by the proximity to its neighbors. In the glassy 
regime, there may be little or no diffusive behavior. Nev- 
ertheless, the approximately constant slope is an indica- 
tor of how much local mobility a given monomer or chain 
can have within its cage. 
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E. The shear modulus from the shear deformation 

We have also calculated the shear modulus directly by 
applying a simple shear deformation to a given sample. 
We start with a cube of side L. For a given strain e^^, an 
affine shear deformation in the x direction is applied (in 
a plane with its normal along the y direction). An atom 
initially in position {x, z) is displaced to {x+Cxyy, y, z). 
The boundaries of the simulation box consequently are 
shifted for Xmin from to Exyy , and for Xmax from L to 
£xyy + L. This is applied as a one-time, instantaneous de- 
formation. The shear modulus is then calculated with an 
off-diagonal element of the stress tensor, once the system 
is equilibrated: 

^ ^ Wxy jexy) - CTxy (0)] ^^y^ 

£xy 

The shearing is applied in the five other directions, 
substituting xy by yz,—xy, etc. Individual samples lack 
symmetry and therefore the various deformations will not 
usually give the same stress components. The evolution 
of the residual stress is monitored, until a stress plateau 
is reached in the deformed system. Deformations of e = 
0.01 to 0.2 were performed. A simple shear deformation 
is more rigorous than a pure shear deformation, which 
relies heavily on the assumption of isotropy throughout 
the system Q. In a pure shear deformation, a sample 
would have been stretched along x by s^y resulting in a 
compression along the two other directions by Sxy/'^ in a 
volume preserving system. 

F. Heat capacity 

Finally, we monitored the changes in the heat capacity 
Cp across Tq- Although other methods exist for calcu- 
lating Cp via MD simulations, the best approach was 
to calculate the potential energy in the NVT ensemble 
(Brownian dynamics) for each temperature. Knowing 
the volume and the pressure at each temperature, we 
can obtain the specific heat per particle Cp in units of 
ks from: 

where H = E + PV is the enthalpy. 

G. The simulations 

Most simulations were carried out on the "Bugaboo" 
Beowulf Cluster at Simon Fraser University, with compu- 
tational times for one complete sample at a given temper- 
ature (ab initio and including all correlations and shear- 
ing) at around 6 days on a single processor. Approxi- 
mately 10 samples of 1050 particles are examined for each 



temperature, with each given sample being used to calcu- 
late approximately 400 "correlators" of length SOOOOdt. 



IV. RESULTS 

A. The glass transition temperature To 

We begin by establishing the glass transition tempera- 
ture Tg, using the common and accessible method (both 
experimentally and by simulation) of examining the vol- 
ume (or equivalently, the packing fraction $) across the 
GT [iO,|43| along the isobar P = 1. The intersection of 
the slopes in the solid and liquid regimes separates two 
regions of nearly constant volume expansivity. It yields 
an unambiguous value of Tq = (0.465 ± 0.005), as shown 
in Fig. ^and Table I. The variation of Cp with respect 
to temperature provides another estimate. As discussed 
in Section lill Fl Cp was calculated from the derivative of 
the enthalpy with respect to temperature along an iso- 
bar (Fig. [SJ. Such an approach has been extensively used 
to calculate Tq and as expected, Cp is characterized by 
a sharp increase around Tq plj . The increase, however, 
takes place over a temperature interval, and it is not clear 
what part of the curve should be used as the location of 
Tq (see Fig. O. The start of the rapid rise occurs at 
0.44 and ends at 0.50. The Tq determined above from 
the change in value of the volume expansivity lies in the 
middle of the rapid rise. The beginning of the rapid rise 
suggests that a major change in the entropy is taking 
place at that point. As we shall see later, this point has 
a special significance. 

B. The viscosity 

Now that we can identify the GT within our sys- 
tem, we can examine the relevant dynamical quantities 
around this point. First, we obtain the expected viscosity 
curve with the characteristic increase over several orders 
of magnitude around Tq, along with the corresponding 
curves of Eqs. (|2Jl and Q obtained using a simple non- 
linear curve-fitting algorithm (Fig. O. We also attempt 
a fit using the MGT approach in Eq. (Q. The results 
are summarized in Table I and follow the same trends to 
those foimd in Refs. 35] and 22] for a similar system. 
The attractive LJ interaction has a shorter cut-off 2.24cr, 
instead of our 2.5cr leading to systematically lower criti- 
cal temperatures. Varnik and Binder have looked at 
the increase in viscosity with lowering of the temperature 
in the same system but with the melt driven by a force 
field. The values of Tmc and Tq arc quite different under 
those conditions, expectedly lower. 

We get an idea of fragility by plotting the log of the 
viscosity rj versus the inverse temperature, a characteris- 
tic Tc-scaled Arrhenius curve for glassy systems (Fig. O. 
Evidence of fragility (deviations from Arrhenius behav- 
ior) becomes apparent as we approach the GT. Typical 
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TABLE I: Relevant temperatures characterizing the glass and 
rigidity transitions: Tq is the glass transition with the asso- 
ciated MCT critical temperature Tmc, the VFT temperature 
To, and the critical temperature Tc obtained from Eq. 
In addition, we present Ti , corresponding to the onset of true 
rigidity, found by extrapolation to Geq in Eq. 11411 . 



To 

0.465 ± 0.005 

Tmc To Tc Ti 

0.51 ± 0.02 0.41 ± 0.02 0.422 ± 0.006 0.44 ±0.01 



experimental data from polymeric glasses generally show 
very fragile behavior from such plots |25^|. 



C. The shear modulus from G{t) 

We have also examined how G{t), the time-dependent 
shear modulus, changes as we cross the glass transition 
(Fig. O. The main difference between solid and liquid 
regimes occurs in the long-time tails of G{t). While a 
KWW fit of Eq. lO can be used for T > Tq, a. power 
law of the form of Eq. H14|) must be applied for T < Tq 
Simply evaluating the value of G'^^ = G{t) at the 

end of our correlation function {t = 150 y^ma^/e) gives 

a good idea of the initial degree of relaxation possible 
in our system, or the shear modulus at short timescales. 
However, we would like to examine the "true," long-time 
relaxation of the glass, which is generally not accessible at 
simulation timescales. As such, we find Geq (at t — > oo) 
using a non- linear fit to the power-law of Eq. H14|l , giving 
non-zero values of Geq beginning at Ti — 0.44±0.01, and 
an approximately linear increase with decreasing temper- 
ature, as shown in Fig. |H1 The decay of G{t) is very slow, 
and most of the information necessary for the fit can be 
extracted from the initial decay region, shown in the inset 
of Fig. 13 Similar power-law behavior has been found in 
colloidal gel simulations We are further convinced 
of this approach by the apparent change of slope of G'^^ 
at the point where Geq becomes non-zero. G'eq begins ac- 
quiring non-zero values at higher T (near Tmc=0.51, the 
idealized MCT temperature). In other words, apparently 
"solid" glasses can have an almost liquid-like response at 
longer times, provided the temperature is just above or 
just below Tq. This is perhaps due to some degree of 
non-ergodicity in the system. We cannot put a precise 
value on the timescale of the long-term rigidity, only that 
it is much greater than that of the simulation. 



D. The shear modulus from a finite deformation 

We also confirm our Geq results by checking whether 
or not they represent the shear modulus /i in its regular 



sense, i.e. the response to a shear deformation, /i, the 
value of G' on our short simulation time scales and the 
value of Geq by fit are shown in Fig. |S1 This preliminary 
comparison is promising, since both Geq and fi{s— 0.1) 
begin acquiring non-zero values at the same temperature 
Ti defined in the previous section: they are both good 
indicators of the onset of true rigidity. Clearly fi{e—0.1) 
is much smaller than Geq ■ while Geq calculates the shear 
modulus via the internal fluctuations, requires exter- 
nal constraints to be imposed, thus altering the system. 
As such, we look for the limit ^(e 0). Fig. IHI shows 
the function fj,(e) for a set of samples at three given tem- 
peratures. This calculation was done by ensuring that 
the stress had reached a plateau before evaluating fj,, via 
Eq. H17(l . Clearly, this glass displays highly non- linear 
behavior and we have, as of yet, no function with which 
to fit the points. The limit e — > is computationally 
prohibitive to reach. The error in /i as seen from Eq. 
rm grows as 1/e. Fig. El does, however, indicate that 
/i(e — > 0) likely tends to Geq. It also suggests that the 
regime of elastic deformation for these glasses is approx- 
imately 2 to 3 %. 



E. The short time diffusion coefficients 

To gain more insight into the structural behavior of the 
melt near the glass transition we look at the short-time 
diffusion coefficients Vm and Vcm of the monomers and 
center of mass of chains respectively, defined in Section 
nil Dl The measurements reflect, on average, how much 
local motion is possible for a given particle. As expected, 
as we decrease the temperature, there is decreased lo- 
cal diffusion as we approach the RT and GT (see Fig. 
2J). In addition we notice that, in the glassy regime, the 
monomers maintain significantly more short-time mobil- 
ity than the chains close to Tq. The chains, on the other 
hand, are more constrained. This is further illustrated by 
the behavior of the ratio T^m /T^cm with temperature. 
Fig. |S1 reveals that monomer motion becomes increas- 
ingly dominant as we approach from above. There 
could be enough free volume for some monomers to move 
short distances . This has been illustrated by the cage 
effect 's^: as the chains become immobilized near Tq, 
the monomers maintain some degree of mobility within 
their "cage" formed by the neighbors. This is, in fact, 
an important precursor to the system becoming not only 
"slow," but also rigid. At Tq the monomers still have 
sufficient mobility to absorb deformations. This may ex- 
plain why the onset of shear rigidity is below Tq. 



V. DISCUSSION 

We have determined the glass transition temperature 
for our system of freely jointed chains interacting with 
van der Waals interactions through the change of slope 
of the packing fraction with temperature, or the change 
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in the value of the volume expansivity. This led to a 
value of Ta = 0.465 ± 0.005. At the glass transition the 
heat capacity at constant pressure Cp exhibits a sharp 
rise with increasing temperature. The above determined 
Tg occurs in the middle of the rise. As seen in Table I 
a number of other temperatures have been obtained in 
our study of the onset of rigidity brought about by the 
kinetic arrest occurring near Tq. The beginning of the 
rise in Cp coincides with Ti = 0.44 ± 0.01 the onset of 
rigidity point as measured by G^q and /i, and the end 
of the rise with Tmc = 0.51 ± 0.02, the MCT critical 
temperature, the divergence point in the viscosity pre- 
dicted by mode coupling theory. These coincidences do 
not appear fortuitous. The rapid rise in Cp suggests a 
temperature interval where the configuration space ex- 
plored by the glass in the making changes rapidly. The 
end of the rise may very well be the temperature at which 
thermally activated processes end and the assumptions 
of idealized MCT become valid. Tmc is expected to be 
higher than Tq 0|. Below Tmc, the liquid can still flow 
with the help of thermally activated hops. The temper- 
ature of rigidity onset was determined by studying the 
behavior of the viscosity r/ and the shear modulus upon 
cooling. The divergence point of rj is obtained by ex- 
trapolation, while the onset of rigidity is measured fairly 
close to the transition. So we will begin with the latter. 
At Tc by all indications the system has no long-time 
shear resistance. Both Geq = limf^oo G{t) and /i, ob- 
tained by an actual deformation of the system, are still 
zero. G'^g = G{t « 150), a short time shear modulus be- 
comes non-zero around Tmc- But Geq, its long time limit 
only sets in at Ti = 0.44 ± 0.01. At the same temper- 
ature n becomes non-zero, although the values of fi are 
not the same as those of Geq for typical deformations of 
e = 0.05 or 0.1 (see Figs. [3ind|Hll- We will elaborate on 
these differences later in the discussion. The important 
point is that they agree on the temperature of rigidity 
onset. On the other hand, the divergence point of the 
viscosity rj was determined using several models, two of 
which incorporate thermally activated processes, which 
become operative below Tmc- They predict divergences 
at To = 0.41 ± 0.02 from a fit to the VFT law (see Eq. 
©) and Tc = 0.422 ±0.006 from the Colby form (see Eq. 
©. One should, however, note that these are obtained 
by extrapolating from the lowest temperature at which 
77 have been measured, 0.49. The true states at Tq or Tc 
are not attained due to the dramatic slowing-down below 
Tc- The extrapolated values are smaller than Ti, but not 
inconsistent with it. They are not expected to be as re- 
liable because they are not measured directly. The fact 
that Ti is smaller than Tc determined from the change 
in volume expansivity raises the question about the true 
location of the glass transition. Does this mean that the 
beginning of the rise in Cp should be used as the true 
onset of the glass transition, or is it, that at Tc, although 
there is structural arrest, sufficient free volume remains 
to allow for small deformations at no cost of energy? 

A number of aspects of this study are worth comment- 



ing on further. First there is always a time scale asso- 
ciated with rigidity, as it is well known that the shear 
response depends on frequency. What we are trying to 
determine is the static shear modulus, or zero frequency 
limit. The long time associated with this modulus, as 
measured, is still short compared with temperatures re- 
lated with aging. We started measuring quantities when 
the system showed no sign of evolving. In other words 
aging is on a different time scale than stress relaxation 
times. An unrelated issue on time scale is the definition 
itself of Tc which is linked to the experimental time scale. 
The glass makers, as mentioned in Section^ use a value 
of 10^'^ Pa-s for the viscosity at Tc- This corresponds 
to the experimental time scale of minutes (or hours) 
while a computer experiment, under the best of condi- 
tions, occurs over nanoseconds, so about 10^^ smaller. It 
is intriguing to note that the VFT fit and the Colby fit 
yield viscosities at Tc of the order of 10^ in our units, 
which with Goo of the order of 10^ gives a relaxation 
time of the order of 10^ at Tc, longer than the time scale 
of our computer experiments, closer to 10"*, or 10^ time 
steps. 

Another point is the nature of the onset of rigidity. 
What happens at Ti? Studies underway li^ show that 
the distribution of displacements on all time scales stud- 
ied is unimodal as previously observed so there is 
no separate long-lived rigid backbone that accounts for 
the resistance to applied shear. Second, the local struc- 
ture of the system remains identical, as expected from 
a GT. In other words, there is no evidence of any clus- 
ters of tightly-packed particles. The average distance to 
the nearest-neighbors of a given particle does not change 
appreciably with temperature. The concept of dynami- 
cal heterogeneity may be more applicable to this system 
[45I li^ , in that the lack of mobile clusters on fairly long 
timescales below Ti accounts for rigidity. Instead of hav- 
ing covalent bonding, the van der Waals rigidity in our 
system may arise from "jamming" constraints produced 
by mechanisms such as the cage effect, which relies pri- 
marily on the presence of stiff chains to prevent motion. 

The fact that the appearance of long-term rigidity at 
Ti ~ 0.44 coincides with the bottom of the characteristic 
"dip" in the heat capacity curve (Fig. is indicative 
of an increase in system stability as the temperature is 
lowered. In network glasses, the number of floppy modes 
in the system determines the jump in Cp (also related to 
the fragihty of the glass) 0. And these floppy modes 
start appearing at Ti. Other approaches have been able 
to associate Cp with the change in entropy of constraint 
breaking |l6l |. Although these approaches rely on physi- 
cal bonds as constraints, the same principles could apply 
to our system. Furthermore, we notice from Fig. |31that 
VcM is very small at Ti. Due to the noise in the data, 
it is not likely that we could resolve the entropic com- 
ponent to the rigidity if it exists. In other words this 
rigidity transition is probably mechanical in nature. 

Finally we would like to comment on the difference 
between Geq and /i, the latter obtained from the applica- 
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tion of an instantaneous simple shear to the system. The 
non-hncar behavior of Fig. I^is further explained in Fig. 
irni Both the initial stress and the degree of relaxation 
are highly dependant on e (Fig. I10|l . For larger shearing, 
the initial stress is smaller and the subsequent decay is 
very fast. A collapse was realized by plotting the stress 
scaled to its initial value versus time scaled with a re- 
laxation time equal to the time to reach half the stress. 
The relaxation times obtained follow the same trend as 
. The shearing is likely irreversible at large strain. In 
effect, large deformations appear to have a substantial 
effect on the structure of the system, perhaps changing 
the shape of the energy landscape and causing the sys- 
tem to be in a different local minimum. Small shear, 
however, does not generally have the same effect and a 
larger residual stress remains as the system tries to return 
to its original configuration (the same energy well). Hav- 
ing described the basic behavior associated with shear- 
ing, we can alternatively examine the local "polymeric" 
contribution to the mechanical relaxation of the glass. 
Considering Fig. ^] we can see that for small e, the 
chains are slightly stretched, but can easily recover their 
equilibrium Rq value, where Rq is the average radius of 
gyration of the chains. For larger strains, the initial de- 
formation of chains is disproportionately high and there 
can be no recovery on intermediate time scales. We can 
interpret the degree of recovery in the chain as the main 
factor characterizing the viscoelastic behavior in poly- 
mer melts (often represented by spring-dashpot models), 
since chain length has a direct, substantial effect on re- 
laxation mechanisms. Furthermore, rigidity should occur 
earlier in a system of longer chains, since each particle 



will have, on average, fewer degrees of freedom. 



VI. CONCLUSION 

Our isothermal compression method has allowed us to 
better examine the state of the melt above and below 
Tg- Our studies of rigidity bring additional insight into 
the nature of the GT. Previous studies on structural is- 
sues related to the glass transition have focussed on the 
divergence of the viscosity. We provide an additional 
perspective by looking also at the emergence of the shear 
resistance. Whereas the point of divergence of the vis- 
cosity is difficult to assess directly because of structural 
arrest, it appears easier to approach the point of onset 
of the shear modulus. That both G^q and /i have the 
same onset is encouraging. The work demonstrates the 
effectiveness of the stress correlators, from which G^q is 
obtained, compared to the application of an external de- 
formation which yields /x. The Tq = 0.465 determined 
from the change in packing fraction lies in the middle of 
the rapid rise of the heat capacity. With decreasing tem- 
perature, there is a short transition period of short-term 
rigidity beginning near Tmc = 0.51 (the end of the rise 
in Cp), presumably due to a lack of ergodicity, followed 
by the appearance of a long-term non-zero shear mod- 
ulus at Ti = 0.44 (the beginning of the rise). I believe 
that we have only begun to explore what shear resistance 
can teach us about the GT: from non-linear mechanical 
properties, frequency dependences, to issues related to 
aging. 
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FIG. 1; Packing fraction <E> as a function of temperature. The 
intersection of the two slopes in the glassy and liquid regimes 
accurately determines To, shown with the arrow and dotted 
hne. 
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FIG. 2: Heat capacity Cp found by numerically deriving the 
enthalpy H with respect to temperature. The dotted line 
shows the location of Ta as determined from Fig. Q Note that 
the drop in Cp as T is lowered ends at T = 0.44, coinciding 
with the appearance of long-term rigidity (see Figl^J. 
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FIG. 3: Viscosity 77 above the glass transition, along with the 
corresponding fits from Eqs. (2) and (3). The dotted vertical 
line marks To obtained previously. 



6.0x10' 










5.0x10'' 




■ 




■ 

■ 


4.0x10' 






1 


1 


3.0x10'' 










2.0x10'' 






■ 

■ 




1.0x10' 
0.0 


- 1 


■ 
• 


■ 

■ • 





0.34 0.36 0.38 0.40 0.42 0.44 0.46 0.48 
T 



FIG. 4: Average short-time diffusion of monomers {T>m) and 
centers of mass of chains {Vcm) around Tq- The substantial 
increase in Vm begins slightly below Tq (indicated by the 
dotted line). 
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FIG. 5; Cage effect as seen by dividing the monomer by the 
chain diffusion constants at short times. The maxiirmm oc- 
curs near To (dotted line), as the cage progressively "closes 
in." 




FIG. 6: Glassy behavior of the polymer melt shown in a typi- 
cal Arrhenius plot, indicating that the glass is primarily frag- 
ile. 
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FIG. 7: Stress autocorrelation function G{t) decay above and 
below Tg, as calculated from Eq. @, along with the corre- 
sponding KWW and power-law fits from Eqs. IjlH^ and 1)14^ . 
Note the very slow power- law decay. Inset: Initial decay of 
the power-law for T = 0.4 shown on a linear plot. 
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FIG. 8: Various approaches to determining the shear modu- 
lus, is simply the value of G{t) a,t t — 150 time units (the 
rigidity on very short timescales), while Geq is obtained by 
fitting the function to Eq. 1141 1 at t — > oo . Finally, /i is cal- 
culated from Eq. 1)17^ via a non-negligible deformation of the 
sample. The inset shows the behavior around the GT (dotted 
line) with the locations of Ti (the onset of long-term rigidity) 
and Tmc (the MOT temperature) indicated by arrows. 
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FIG. 9: Shear modulus /i as a function of shear for three 
temperatures below the GT. We interpret the Geq as the zero- 
shear limit to check whether these two methods agree. Note 
how the statistics get much worse at small e. 




FIG. 10: Decay of of the shear stress axy after various initial 
shear deformations at T = 0.3. We can see both the initial 
stress and the shape of the subsequent decay in the system. 
Note that for large e there is very little initial stress imparted, 
followed by a fast decay. 
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FIG. 11: Evolution of the average radius of gyration. Eg, of 
the polymer chains, after a largo (e = 0.2) and sniall (r = 
0.05) deformation at t = together with the undeformed, 
"equilibrium" Rq (e = 0). 



